Can we avoid dark energy? 
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The idea that we hve near the centre of a large, nonlinear void has attracted attention recently 
as an alternative to dark energy or modified gravity. We show that an appropriate void profile can 
fit both the latest cosmic microwave background and supernova data. Ifowever, this requires either 
a fine-tuned primordial spectrum or a Hubble rate so low as to rule these models out. We also show 
that measurements of the radial baryon acoustic scale can provide very strong constraints. Our 
results present a serious challenge to void models of acceleration. 

PACS numbers: 98.80.Es, 95.36.+X, 98.65.Dx 



Introduction. — The last decade has seen the solidifi- 
cation of the Standard Model of Cosmology (SMC; see, 
e.g., [1]), which contains about 75% dark energy driv- 
ing the acceleration of a flat, homogeneous and isotropic 
Fricdmann-Lcmaitrc-Robcrtson- Walker (FLRW) back- 
ground. A very broad range of observations arc con- 
sistent with the SMC, including cosmic microwave back- 
ground (CMB) (sec, e.g., [2]), Type la supernovac (SNc 
la) (e.g. [3]), baryon acoustic oscillations (BAO) (e.g. [4]), 
integrated Sachs-Wolfe (ISW) effect correlations (e.g. 
[5]), weak lensing studies (e.g. [6]), etc. 

Given this impressive convergence of observations, it 
would be very surprising if the data still fit an alternative 
model, which lacked two of the main planks of the SMC. 
Yet just such a radical proposal has attracted consider- 
able attention recently. The idea is to drop the dark en- 
ergy and the Copernican principle, and instead suppose 
that we are near the centre of a large, nonlinearly under- 
dense, nearly spherical void surrounded by a flat, matter 
dominated Einstein-de Sitter (EdS) spacetime ([7]; see [8] 
for a brief review). By tuning the radial void profile, it 
is possible to match the luminosity distance-redshift re- 
lation of the SMC [9] , so if these models are to be ruled 
out, we need more than just the SN data. Recently, con- 
straints from the CMB acoustic scale, BAO scale, and 
Hubble rate have been placed on void models, although, 
remarkably, they have not yet been ruled out [10, 11]. 

We must stress that void models contain two elements 
which appear extremely unlikely within the standard cos- 
mological framework, and for which no viable explana- 
tions have been proposed. Voids of the size required to fit 
the SN data (hundreds of Mpc to Gpc scales) correspond 
to fluctuations of very many a in standard structure for- 
mation scenarios [12]. Also, we must be very close to the 
void centre to avoid a large CMB dipole [13]. 

However, as overwhelming as these difficulties may ap- 
pear, they are essentially philosophical in nature, and so 
it would be important to rule these models out on the ba- 
sis of observations. With this goal in mind we confront 
void models with several sets of current data, providing 
three main advances over previous studies. Firstly, we al- 
low for a wide range of void profiles, employing a spline 



parameterization. Secondly, we carefully calculate the 
CMB anisotropy spectrum that a void observer would 
see and confront the full spectrum with recent data sets. 
Finally, we show that the radial BAO scale is a powerful 
discriminator between void and standard models. 

Specifying the void model. — We model the contents of 
the Universe at late times as pressureless matter, with 
density p. Observations arc consistent with isotropy, so 
we place the observer at the centre of spherical symmetry. 
The exact solution to Einstein's equations in this case is 
known as the Lemaitre-Tolman-Bondi (LTB) spacetime 
[14]. It is described by two free radial functions, which 
correspond to the growing and decaying modes in the 
limit of small perturbations about FLRW [15, 16]. As 
stressed in [16], it is crucial to consider only voids with 
vanishing decaying mode (i.e. uniform "bang time"), if 
we are to be able to specify the initial conditions (ICs) 
for perturbations at early times. This is because in this 
case the void itself will be a small perturbation from 
FLRW at early times, and standard inflationary ICs can 
be specifled on top of the void. Since our analysis will 
include the BAO perturbation scale (which is set before 
recombination) evaluated inside the void, we must as- 
sume vanishing decaying mode in our work. Thus only 
one radial function is required to specify the void model. 

The LTB spacetime can be described by the metric 
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Here Y{t, r) and K{r) are determined by the exact LTB 
solution (see, e.g., [16]) once the single free radial func- 
tion is specified, and Y' = dY/dr, where r is a comoving 
radial coordinate. Errors due to ignoring radiation are 
inevitable with the LTB solution, but wc estimate our 
results are accurate to the percent level or better [17]. 

We chose to define the radial profile in terms of the 
comoving density perturbation, 8p{ti,r) / p{ti), specified 
at the early time t;. To be confident that no impor- 
tant regions of "profile space" are missed, we fit a three- 
point cubic spline to the initial density fluctuation 5pj = 
Sp{ti,rj), where j = 1, 2, 3. We fix ri =0 and enforce the 
void to be smooth at the origin and to smoothly match 
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to EdS at large r by setting Sp'{ri) = Sp'{r3) = 5pz = 0. 
This leaves a total of four free parameters: the density at 
the origin, 5pi, the density and radius at the midpoint, 
5p2^ r2, and the radius at which we match to EdS. We 
have checked that additional spline points do not signif- 
icantly improve the fit to current data. We consider two 
classes of profiles: "constrained," for which we impose 
/ 5p{ti)r'^dr < 0, and "unconstrained," which are free. 

A void model is completely specified by the radial pro- 
file, the Hubble rate at the void centre today, Hq, the 
radiation density today, which is fixed by the CMB mean 
temperature, Tq = 2.725 K [18], and the baryon fraction 
fh = pb/ Prn- Outside the void we asymptote to EdS. 

In Fig. 1, wc plot several void profiles on the observer's 
past light cone, in terms of the local density parameter 
= 2AitGp/9^ , where 9 is the comoving expansion [16] 
(this definition reduces to the standard density parameter 
in the FLRW case). The profiles are sampled from a 
Markov-Chain-Monte-Carlo (MCMC) process, fitting to 
SN-f CMB data as described below. In these plots, values 
< 1 correspond to the central void, while il^^ > 
1 indicates an overdense shell region. The constrained 
voids tend to be essentially "compensated," while the 
unconstrained tend to be strongly "overcompensatcd." 
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FIG. 1: Local density parameter Q^^'^ versus redshift z for 
constrained (left) and unconstrained (right) voids. The pro- 
files are sampled from MCMC chains and the grayscale level 
indicates the relative likelihood. 

Cosmic microwave background. — The CMB anisotropy 
power spectrum, Cg, contains much information, and so 
can potentially provide strong constraints on void mod- 
els. Two main factors determine the detailed shape of 
the Ci spectrum: the primordial perturbation spectrum, 
which, in the simplest models of inflation, is close to scale 
invariant and essentially featureless; and the local physics 
during recombination, which is determined by the com- 
position of the Universe at that time. The angular scale 
at which features of fixed physical scale appear is deter- 
mined by the physical circumference of the observer's last 
scattering surface (LSS), which is affected by the geom- 
etry of the Universe and the time of observation [19]. 

Once the void model has been specified, it is straight- 
forward to calculate the Ce spectrum at all but the largest 
scales. The anisotropics generated at the LSS cannot be 
affected by the void for an on-centre observer, because 



of the isotropy of the LTB background [20]. Our proce- 
dure to calculate the 's is simply to find the parameters 
for an effective EdS model with the same physics at re- 
combination and LSS physical circumference as the void 
model. Then those effective parameters can be fed into 
public CMB anisotropy codes — we used CAMB [21]. 

Explicitly, we numerically integrate along the past 
light cone from the void centre today {r = z = 0) 
to a redshift Zm to find the corresponding coordinates 
itm,rm) = (t(zm),r{zm)) using the exact LTB relations 

dr _ VI- K dt _ -Y' 

dz~(i + z)Y'' dz^{i + z)Y'' 

We choose Zm large enough so that is sufficiently far 
outside the void that spatial curvature (and shear) is 
negligible there, but not so large that radiation is im- 
portant at tm at background level. In practice, values 
Zm — 100 meet these criteria. We also evaluate the 
Hubble rate, Hm, at z„j using the exact LTB solution. 
Finally, we set the spatial curvature to precisely zero, 
and integrate back up the light cone into EdS to co- 
moving coordinate r^'^^ = 0, using the FLRW relation 
dz^'^^/dr^'^^ H^'^^{z). The correct LSS circumference 
is ensured by using the relation Y{tm,rm) = «(''m'^^)^m'^^j 
for FLRW scale factor a, to set the IC for this integration. 
The result of integration, z^^, allows us to calculate the 
effective EdS mean temperature and Hubble rate via 

/-i I ^EdS\ 3/2 

T^^^ = To(l+zr), Hi^' = Hm (-7^) ■ (3) 

The parameters T^'^^ , H^<^^, n,n = I, = ^ 0, 
and fb are then fed into CAMB to calculate the C( spec- 
trum which would be observed in the specified void. Note 
that the effective parameters Tq'^^ and Hq*^^ will gener- 
ally differ from the true temperature and Hubble rate 
observed at the void centre today, Tq = 2.725 K and Hq. 

Importantly, we find that the effective set Tq*^^ = 
3.372 K, ^ 5i.0kms-iMpc~i, and = 0.165 pro- 

duce a Ci spectrum that matches that of the Wilkinson 
Microwave Anisotropy Probe (WMAP) best-fit A model 
[2] at all but the largest scales. A void model must have 
effective parameters close to these if it is to fit the CMB. 

Data fitting and results.— We used COSMOMC [22] to 
create MCMC chains to estimate confidence limits on 
the parameters. Along with the four void spline quan- 
tities, the basic cosmological parameters used in the fit 
arc fb, Hq, and the amplitude As and spectral index Us 
of adiabatic initial perturbations [20]. 

SNe la provide important evidence for acceleration 
within the standard FLRW framework, and hence we 
must ensure that our voids fit these observations. We 
used the recent Union compilation [3], consisting of 307 
SNe with z = 0.015-1.55 [23]. We also used CMB data, 
namely the 5-year WMAP results [2] and four exper- 
iments at higher resolution: ACBAR [24], Boomerang 
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FIG. 2: Likelihood contours at 68% and 95% confidence for 
constrained (blue, solid contours) and unconstrained (red, 
dashed contours) voids. 
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FIG. 3: Residual Ce's (in fjK ) from the best-fit A model for 
constrained (top) and unconstrained (bottom) voids. Binned 
WMAP data are shown by the red (thick) error bars, and 
ACBAR data by the blue (thin) bars. 



[25], CBI [26], and QUaD [27]. We also applied a conser- 
vative prior that il}^^ > 0.1 at the void centre. 

In Fig. 2 wc show a selection of 2D likelihood sur- 
faces for the constrained and unconstrained void mod- 
els. The parameters are very different in each case. 
For constrained voids, the effective temperature, Tq*^^ = 
2.760 ±0.008 K, is similar to Tq. This T^'^^ is far too low 
to provide a good fit to the observed CMB spectrum — 
we find Ax^ = 162 between the constrained void and A 
model for CMB-I-SN data, with almost all of this differ- 
ence coming from the CMB. This poor fit also leads to 
very low fb = 0.100 ± 0.001 and n, = 0.88 ± 0.01. 

For the unconstrained voids, however, much higher ef- 
fective temperatures are possible, due to the geometrical 
effect of the overdense shell. These temperatures are suf- 
ficiently high to fit the CMB well— the fit to CMB+SN is 
actually slightly better than A, with Ax^ = —1.4. How- 
ever, this requires an unusually low local Hubble rate of 
iJo = 44±2 kms~^ Mpc~^, as Fig. 2 shows. Recent local 
estimates range between 57 and 79 kms~^ Mpc~^ at la 
[28] , and so this class of void model is ruled out at high 
confidence. 

In Fig. 3 we present the residuals of the void C^'s from 
the best-fit A model. The best-fit void model is shown 
by the black curve, along with 100 other spectra sampled 
randomly from the MCMC chains. The grayscale level 
indicates the relative likelihoods. It is clear that in the 
constrained case one cannot fit the CMB data without 
introducing fine-tuned features to the primordial spec- 
trum, since the physics at recombination is wrong. 

Baryon acoustic scale. — The physics before recombi- 
nation imprints a fixed comoving scale into the matter 
power spectrum in the form of the BAO scale, carrying 
much useful information about the geometry and expan- 
sion history assuming an FLRW background (see, e.g., 
[29]). It is therefore important to assess the usefulness of 



BAO data in constraining void models for acceleration. 

The first step is to evaluate the physical BAO scale, 
if^'^ (also called the sound horizon at the drag epoch), 
at some time ti early enough that the void background is 
well approximated by FLRW. To do this, we find an ef- 
fective EdS model with the same physics at early times as 
the specified void model using the same procedure used 
to calculate the void C^'s, except that it is not necessary 
here to match the effective and true LSS circumferences. 
Then wc calculate if^'^ using the fitting function from 
[30] applied to the effective EdS parameters. Next, if^'^ 
is propagated up to redshift z on the void observer's past 
light cone using the LTB metric. The background shear 
causes the physical BAO scales in the transverse and ra- 
dial directions to differ; they are given, respectively, by 



BAO 
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(4) 



Here Y{z) = Y{t{z),r{z)), and similarly for Y'{z). BAO 
observations are often expressed as a BAO length scale 
today, but such values necessarily depend on the assumed 
background. Instead, a model-independent expression of 
the transverse and radial BAO (RBAO) scales is in terms 
of the corresponding angular and redshift increments. 
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Reference [31] emphasized the importance of distinguish- 
ing radial and angular scales in this context. 

Importantly, Az^^*-* contains two factors that rein- 
force each other in the peripheral void region. The quan- 
tity Y'{z)/Y'{z) is the expansion rate in the radial di- 
rection, which is suppressed in the overdense periph- 
ery. This, in turn, results in a suppressed RBAO scale, 
l^^'^{z), in the periphery. The net effect is a strong sup- 
pression 

of Az^AO in 

this region compared with the stan- 
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dard A illustrated in Fig. 4 for a set of profiles 

from the same MCMC chain as we obtained above [32] . 




the Copernican principle. 
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FIG. 4: RBAO scale Az^*° for constrained (left) and uncon- 
strained (right) voids. Also shown are the RBAO data from 
[33] (error bars), best-fit A model (blue, dashed), and best-fit 
void including RBAO constraints (green, dot-dashed). 



Recently, RBAO measurements have been presented 
[33] for the observable Az^^*-* rather than the model- 
dependent Z^^*-*. Figure 4 shows that the new data are 
in excellent agreement with standard A, but strongly at 
odds with the likeliest void models. Including the new 
RBAO data points worsens our best fits by Ax^ = 16 
(constrained profiles) and 9.8 (unconstrained). Such 
large A^^ values for only two extra data points show 
that the new RBAO technique is already a considerable 
obstacle to void models. 

Discussion. — We have concentrated here on con- 
straints from SN, CMB, and RBAO data. In future work 
[17] we will also consider the independent constraints 
from the amplitude of matter fluctuations [16]. These 
techniques all rely on observations confined to our past 
light cone. Very promising are measurements that sam- 
ple the interior of the light cone, via spectral distortions 
[34] or the kinematic Sunyaev-Zel'dovich effect [35] (see 
also [36] for a distinct approach). 

Nevertheless, we have already shown here that an ap- 
propriate overcompensated void profile can fit both the 
SN and CMB data, but at the expense of an Hq so low 
that it can be ruled out. Constrained voids provide a 
better Hq, but a very poor fit to the CMB without fine 
tuning of the primordial spectrum. We stress the signif- 
icance of this result: it is an extraordinary achievement 
of the SMC that it predicts the detailed shape of the 
Ce spectrum using only a few parameters. Losing this 
predictive power and requiring a fine-tuned primordial 
spectrum is a severe price to pay for the allure of A = 0. 

We also showed that early RBAO results already im- 
pose very strong constraints on void models. In partic- 
ular, RBAO poses a serious challenge to the constrained 
profiles which is free of any subjectivity that some may 
argue is inherent in the above fine-tuning argument. Fu- 
ture BAO surveys such as PAU [37] and BOSS [38] will 
provide improved precision out to greater redshifts, plac- 
ing unprecedented constraints on inhomogeneity. In this 
era of precision cosmology wc can in fact begin to test 
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